Theory of the non-equilibrium quasiparticle distribution induced by Kondo defects 



O 

o 

(N 

U 

Oh- 
< 
On 



J. Kroha a and A. Zawadowski h ' c ' d 
a Institut fur Theorie der Kondensierten Materie, University of Karlsruhe, POB 6980, D-76128 Karlsruhe, Germany 
b Department of Physics and c Research Group of the Hungarian Academy of Sciences, Technical University of Budapest 
d Solid State and Optical Research Institute of the Hungarian Academy of Sciences, POB 49, H-1525 Budapest, Hungary 

(February 1, 2008) 

It is shown that in resistive nanowires out of equilibrium containing either single- or two-channel 
Kondo impurities the distribution function f(E, U) obeys scaling behavior in terms of the quasipar- 
ticle energy E and the bias voltage U. The numerically calculated f(E, U) curves explain quantita- 
tively recent experiments on Cu and Au nanowires. The systematics of the impurity concentration 
Cimp extracted from the comparison between theory and results on various Cu and Au samples 
strongly suggests that in these systems the scaling arises from magnetic Kondo impurities. 
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Electronic interactions in solids are usually investi- 
gated by means of linear response or by spectroscopic 
measurements which essentially probe the system in ther- 
modynamic equilibrium. In an important series of exper- 
iments 0,| the Saclay group demonstrated that unique 
information about the energy dependence of the inter- 
actions in a mesoscopic wire can be extracted from the 
shape of the distribution function f x (_E, U) of quasiparti- 
cles (qp) with energy E at a position x in the wire, when 
a controlled non-equilibrium situation is established by 
applying a finite transport voltage U . In resistive Cu 
and Au nanowires out of equilibrium the theoretically 
expected double-step form of f x (E,U) || was found 
to be rounded such that it obeys the scaling property 
f x (E,U) = f x (E/eU), when U exceeds a certain energy 
scale By a phenomenological analysis within 2nd 

order perturbation theory, the origin of the scaling be- 
havior was traced back to an effective electron-electron 
interaction v(u>) which scales with the energy transfer uj 
as v (uj) oc 1/uj §. Such an lu dependence implies, in 
particular, that v(u>) has no essential momentum depen- 
dence and, hence, should be of local origin. Moreover, 
within the perturbative treatment it implies a logarith- 
mic divergence of the energy relaxation rate at the Fermi 
energy Ep. The latter has generated substantial interest 
because of the possible relation to the problem of dephas- 
ing saturation Q in mesoscopic systems. 

Anomalous low-energy behavior of local origin can be 
induced by the Fermi surface singularities characteris- 
tic for Kondo type systems Inelastic scattering by 
Kondo impurities was discussed in fj|. Based on these 
considerations, the single-channel Kondo (ICK) Q and 
the two-channel Kondo (2CK) effect 0, possibly pro- 
duced by degenerate dynamical defects ||, have been 
proposed as the origin of the anomalous energy relax- 
ation. In this Letter we show that a very small concen- 
tration Ci mp of either ICK or 2CK impurities leads to the 
observed scaling behavior of f x (E, U), when ell exceeds 
an intrinsic energy scale eU* which is essentially equal to 



the Kondo temperature Tk- The numerical results are in 
excellent quantitative agreement with the experimental 
curves 0j2|, with Ci mp the only adjustable parameter of 



the theory. A detailed analysis suggests that the scaling 
behavior in Cu and Au wires is due to magnetic Kondo 
impurities. 

Let us first set up the general formalism for calculat- 
ing f x (E, U) in a resistive nanowire of length L, subject 
to the boundary conditions that the left (x = 0) and 
the right [x = L) leads are in equilibrium at their re- 
spective chemical potentials, i.e. f x =o(E,U) = f°(E), 
fx= L (E, U) = f°{E + eU), where f°(E) = l/(e E ' T + 1) 
is the Fermi distribution (ks — 1). The lesser (<) 
and the greater (>) conduction electron Keldysh Green's 
functions read G x (p,E) — —2irif. x (p)ImG x (p,E) and 
G>(p,E) = 2m[l - f x (p)]l^G r x (p,E), respectively, 
where E and p denote energy and momentum of the 
quasiparticles in a small volume centered around x, in 
which the external fields may be considered constant. A 
superscript r indicates here and in the following a re- 
tarded propagator. In a disordered electron system with 
diffusion constant D the stationary quantum Boltzmann 
equation for the distribution as function of the qp energy 
E takes the diffusive form 01, 



DV 2 J x (E,U)=C{f x (E,U)} 



(1) 



The collision integral C is expressed in terms of the self- 
energies for scattering into (<) and out of (>) states 
with given energy E (N Q — density of states per spin) as 

' ^p<^)G>(p,£)-S>(E)G<(p,i;)]. (2) 
p 
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In the absence of any interactions (C = in Eq. (0)) the 
distribution function has the double-step shape, 



fx(E,U) 



X f°(E + eU)+(l-^)nE) . (3) 



L ' \ Lj 

For a small concentration of Kondo defects Ci mp , in ad- 
dition to the static impurities, the conduction electron 
selfenergy is given in terms of the single-particle t-matrix 
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FIG. 1. a) Diagrammatic representation of the NCA equa- 
tions Solid, dashed and wiggly lines denote the con- 
duction electron, the local (pseudo)spin and the auxiliary bo- 
son propagators, respectively, b) Leading contribution to the 
inelastic spin relaxation rate induced by a finite bias voltage 
U (see text), c) Decomposition of the NCA single-particle 
t-matrix t to show that it includes the Kondo induced elec- 
tron-electron vertex of 0(J 3 ) and 0(J 4 ). 

of the defect, t$(E), as s| = c imp t^. t^{E) mediates 
energy transfer between electrons in that it couples the 
dynamical defect both to the in or outgoing electron and 
to intermediate particle-hole pairs. The elastic scattering 
parts of tx(E) cancel each other exactly in C. We em- 
phasize that, apart from the assumption of small Ci rnp , 
the present formulation, Eqs. (0), (|2|), contains no ap- 
proximations, once the i-matrix is known. 

As pointed out in Ref. |J, the precise energy de- 
pendence of the electron-electron vertex is crucial for 
whether or not f x (E, U) obeys a non-equilibrium scaling 
property, but has been notoriously difficult to calculate 
for the Kondo problem. While 4th order (unrenormal- 
ized) perturbation theory yields the correct power law 
for scaling [^J, partial summations of logarithmic terms 
give corrections violating scaling H|l(J. However, such 
summations are valid only for T, E 3> Tr Chapt. 3, 
while the non-equilibrium situation (elJ > Tk, T <C Tk) 
may involve all energies T < E < eU . 

Therefore, we use the slave boson formalism, where 
certain exact properties of the auxiliary particle propaga- 
tors are known J]~l|— |T^] . To describe 1CK as well as 2CK 
impurities we use the SU(N)xSU(M) Anderson impurity 
model in the Kondo limit, denoting the spin degeneracy 
by N and the number of identical, conserved conduction 
electron channels by M. Following the notation of Ref. 
lUll, the hamiltonian reads, 

H = H + e d flfa + (f$ b ™ c pm* + h.c.) , (4) 

cr p,m,(T 

subject to the operator constraint Q — 2 ff /]/i7 + 

Em&flAn = I- H = Ep,rn, ? £ P C l™a C pm<y describes 

the conduction band. The auxiliary fermion and boson 
operators, /J, b\- n , create the local defect in its quan- 
tum state a or in the unoccupied state, respectively. In 
the 1CK case (N = 2, M = 1) of a magnetic Ander- 
son impurity a denotes spin, and m = 1 has no rel- 
evance. For a 2CK defect (N = 2, M = 2), a is 



identified with a pseudospin, e.g. the parity of the lo- 
cal defect wave function, and m = 1,2 is the conduc- 
tion electron spin, acting as the channel degree of free- 
dom. The equilibrium Kondo temperature of the model is 
~ E F {NN J)W N )e-^ NN ° J \ with J = \V\ 2 /\e d \ 
the effective spin exchange coupling. The bare auxil- 
iary particle propagators read G^ (w) = l/(o; + iO) 

and G r b ^°\ijS) = l/(w + Ed + iO). Here we have gauged 
the zero of the slave particle energy such that the pole of 
G r j is at uj = 0. 

The numerical evaluations of physical quantities will 
be done within the non-crossing approximation (NCA) 
which is shown diagrammatically in Fig. a). The corre- 
sponding equations for the auxiliary fermion and boson 
selfenergies £y , E^- in non-equilibrium read |]l4| , 

where T — ttN \V\ 2 is the effective hybridization, and 
Gjr (e) = ^2 p G$ (p,e). This set of selfconsistent, non- 
linear equations is closed by the Kramers-Kroenig rela- 
tions, Gj b {uo) = — J de/(2m) G^ b (co)/(oj — e + iO), which 
follow from causality and the fact that the auxiliary par- 
ticle Green's functions have only forward in time propa- 
gating parts. Within NCA the single-electron i-matrix 
due to the Kondo impurity is 

By writing the renormalized auxiliary propagators as 
the bare ones with selfenergy insertions and using 
|V| 2 Gr (w rs 0) = J, it is seen that the present formula- 
tion includes the Kondo induced electron-electron vertex 
of 0(J 3 ) and 0(J 4 ), with the points 3 and 4 connected 
(Fig. |c)). The exchange diagrams (points 3 and 2 con- 
nected) are not included. However, this does not change 
the scaling properties (see below). By self-consistency, 
NCA goes beyond this two-particle scattering approxi- 
mation considered in Ref. ||. 

The NCA in non-equilibrium includes both an inelas- 
tic spin relaxation rate and scaling behavior in terms of 
the applied bias. It is instructive to investigate these 
properties for a single Kondo impurity before we present 
the numerical solutions. The Kondo scale Tk is influ- 
enced by the step heights in f x (E, U), Eq. (|J), and the 
size of the corresponding logarithmic terms in the Kondo 
vertex, which is reduced compared to equilibrium. This 
leads to a suppression of Tk, e.g. in the middle of the 
wire (x/L = 1/2), 

I — — e[/»T (0) t(°) 2 

T K = V(ef//2) 2 + 2 - eU/2 ~ * . (8) 
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At an arbitrary position x/L, for eU we have 



T K = T { ° } l,v /{eU)W- x , where r\ = max[x/L,l - 
x/L]. At the same time, there is an inelastic spin relax- 
ation rate l/r s , since in the non-equilibrium electron sea 
(Eq. (||)) there is finite phase space available for scat- 
tering even at T — 0. Technically, this relaxation rate 
appears as the imaginary part of the pseudofermion self- 
energy, Yj r j(Lo — 0), which carries the local spin degree of 
freedom. To leading order in J it is obtained by inserting 

>(0) 

the bare propagators G^ b 



in the diagram Fig. |l| b), 
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2ttMN-(i - - )(N J) 2 eU . (9) 

This is analogous to the well-known Korringa spin relax- 
ation rate §], with T replaced by eU. Solving Eqs. (|), 
(0) selfconsistently in the complete range of validity of 
NCA, Tk < eU <C Ep, we find that beyond perturba- 
tion theory 1 /r s depends on eU and Tk only, 

\H MtN (?p)eU , (10) 



1 x / x 
~~ IV _ L 



where the universal function Hm,n(v) — > 7rM/[2iVln (y)] 
for y ^> 1 |t5), in accordance with Ref. S. Inserting 
l/r s into the pseudoparticle propagators, it cuts off all 
logarithmic contributions of perturbation theory. Thus, 
the low-T scale of the non-equilibrium Kondo system 
is T = max[TR-(eJ7), l/2r s (e/7)]. The crossover from 
the Kondo (Eq. (|8| )) or T limited life time to the inelas- 
tic scale (Eq. (|10|) ) occurs as function of eU at a bias 
eU*. It follows from the universality of l/r s that ei7* 
is only a function of T^\ i.e., for dimensional reasons, 
eU* = A MrN T$\ Numerically we find A 1>2 = 1.48±0.08 
and A 2 . 2 = 1.39 ± 0.05. For eU > eU* T K has lost its 
relevance, and for eU > Well* one has with good accu- 
racy X/t s oc eU, when eU is varied by a factor of ~ 4, 
wherein the ei7 dependence of the log terms is weak. 

To investigate scaling of f x (E, U) we must consider the 
energy dependence of the exact pseudoparticle Green's 
functions, G^ b (u>), from which all other physical quanti- 
ties are derived. It is known that in equilibrium at T = 
it is determined by an infinite logarithmic series which 
results in power law behavior, G"j(uj) oc iQ(±uj)\uj\~ af , 

Gf(uj) oc ie(±ui)\cj\- ab for uj < T { ° ] . The expo- 
nents a/, «b are due to an orthogonality catastropy in 
the auxiliary propagators and have characteristic values 
a f = a b = 1/2 for the ICK and a f = M/(M + N), 
a b = N/(M + N) for the 2CK fixed point of the model 
Eq. (||) (Kondo limit) We can exploit this knowl- 

edge to determine the frequency dependence of G"j b {oj) 
away from equilibrium without explicitly summing up 
the logarithmic series. At finite bias eU 3> Tk this series 
consists of similar terms as in equilibrium, however with 
three modifications: (i) Because of the inelastic relax- 
ation rate l/r s = 27 all frequency arguments are shifted, 
u) = uj + ij. (ii) Gf(oj) has a singularity at to — + 17, 



but there are two singularities in Gf (uj) at to = + ib"f 
and at uj — eU + itry, where 6 is a numerical factor, (in) 
Each frequency integral involving Gf(e + uj), like, e.g., 
in Eq. (|), carries a prefactor M, and each of the two 
singularities in Gf gives a singular contribution of equal 
weight at the external frequency u> = 0. This can be seen 
as an effective doubling of M . Points (i)-(iii) can be 
verified by iterating Eqs. (||), (||), starting from the bare 
propagators Gy b . As a result, we obtain at x/L = 1/2 
damped power law behavior for the auxiliary propagators 
in non equilibrium, 



Gj (uj) oc ilm 
G^ (us) oc ilm 



l-f 1/2 (E, U) 
(uj + ij) a f 

l-f 1/2 (E,U) ^ 

(uj + Hry) a 'i ^ (uj-eU + ilry) a '» 



(11) 



l-fi/ 2 (E,U) 



The exact exponents [ [l3[ in the non-equilibrium situ- 
ation (with M — > 2M in the logarithmic series) are 
a' f = 2M/(2M + N), a' b = N/(2M + N). The uj de- 
pendence Eq. ([n]) should extend from uj = up to the 
smallest energy scale of the model, i.e. for eU > eU* up 
to uj — eU, since in this case the Kondo scale has dis- 
appeared. The behavior described above is confirmed by 
our numerical NCA solutions. For x/L^0oix/L^ 1 
the solution crosses over to the equilibrium one, as ex- 
pected. The modification of the exponents a'^, a' b com- 
pared to their equilibrium values is reminiscent of a dou- 
bling of the channel number due to the two Fermi edges. 
It remains to be seen whether a strong coupling region 
(T^ 0) < eU < eU*) can be realized where such behav- 
ior can be observed in the presence of l/r s ~ 0(eU). 
The latter was neglected in Ref. |l(|. Here we are in- 
terested in scaling at large bias (eU 3> elf*). Inserting 
the power law forms Eq. (0) into Eqs. (|)-(0), divid- 
ing Eq. (|) by (eU) a 'f and Eq. (|) by (eU)< , and using 
the exact result a'j + a' b = 1, it is seen that the NCA 
equations contain only dimcnsionless energies, e/ eU etc. 
Power counting arguments Jk| show that this is repro- 
duced in arbitrary selfconsistent order in T beyond NCA. 
In the presence of a finite concentration Ci mp , f x (E, U) is 
determined by the selfconsistent coupled set of equations 
(0), (|) and (|)-(0). It follows that the solution obeys 
scaling, f x (E,U) = f x (E/eU) for eU > eU*. Our nu- 
merical solutions show scaling within a factor of 4 to 9 in 
eU, depending on parameters, wherein log corrections to 
1/t s oc eU, Eq. (H]), are small. Note that the power law 
behavior Eq. (|ll ) and the fact that the low-energy cutoff 
1/t s itself is proportional to eU (up to small log correc- 
tions) cooperate to produce scaling. For ell < 10eC/* we 
find deviations from scaling, because then the latter con- 
dition is no longer fulfilled. This provides for T <§; Tj^' 
a rough estimate, and for T > T^ an upper bound on 
T { ° ] - in the experiments @,| T < T^ <C eU. 

For the numerical evaluations we assume magnetic 
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E/eU E/eU E/eU 

FIG. 2. Non-equilibrium distribution functions for various Cu and Au samples. Black lines: experimental results; Cu: [1], 

Au: [2]. Open circles: theory for eU 3> ell*. Deviations from scaling at smaller ell [1] are also reproduced by the theory (not 

shown). Fitted d m p values are indicated. The insets show the difference between the experimental and the theoretical curves. 



(1CK) impurities (2CK impurities give very similar re- 



sults) and take T^ 0) w 0.1 K in Cu and T^' w 0.5 K in Au 
wires (corresponding to N a J — 0.041 and N a J = 0.048, 
respectively), consistent with the above estimate and 
with independent estimates of for these samples [|| . 



(0) 



After t£ } is fixed, 



is the only adjustable parame- 
ter of the theory. The results for f x (E, U), as measured 
by a tunnel junction attached to the wire, are shown in 
Fig. ^. Excellent quantitative agreement with experi- 
ments |l],^ is obtained for all samples. In Au wires the 
fitted values of Ci mp are consistent with (although some- 
what higher than) independent estimates of the magnetic 
impurity concentration [pi, considering the roughness of 
both estimates. This suggets that the scaling behavior of 
f x (E, U) in the Au samples is due to magnetic (1CK) im- 
purities. Furthermore, in all Cu samples the fitted Cj mp 
is ~ 10 2 times smaller than in Au. This systematics is 
in accordance with Ci mp estimated from the plateau in 
the T dependence of the dephasing time t v in similarly 
prepared samples |2p7|]. 

In conclusion we have shown that single- or two- 
channel Kondo impurities in quantum nanowires induce 
scaling behavior of the non-equilibrium distribution func- 
tion f x (E,U) at a bias ell exceeding an energy scale 
eU* re T { °\ The results give a detailed explanation of 
related experiments. In the small bias or strong coupling 
regime (Tk > eU), 1CK and 2CK impurities must show 
qualitatively different behavior, as the former become po- 
tential scatterers with frozen spin dynamics, contrary to 
the latter with (ideally) non-zero entropy at T = 0. The 
quantitative comparison between the present theory and 
experiments suggests that in Au and at least partially in 
Cu nanowires both the scaling of f x (E, U) [j]j|] and the 
plateau in the low-T dephasing time t v |17] are due to 



magnetic Kondo impurities. A unique test for magnetic 
impurities will be measuring f x (E, U) in a magnetic field. 
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